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Abstract. The first theoretical study of a dimers adsorption process at homogeneous 
surface is presented. By using the RSA algorithm, we show example monolayers, 
discuss estimations of random jamming coverage and measure the surface blocking 
function, which could be used for calculating real systems kinetics. We also found the 
correlation function for generated coverages and analysed orientational ordering inside 
the adsorbed monolayer. Results were compared with theoretical and experimental 
data. 
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1. Introduction 

Adsorption (in particular irreversible adsorption) of bio-particles at interfaces plays 
an extremely important role in biotechnology, medicine, chemistry and environmental 
technology. Some common examples are: paper production [lj, particle deposition 
controlling [2] and selective deposition of cells and viruses (31 H] . Simultaneously, there 
has been much interest focused on patterned surfaces with regular shape features like 
circles and dots, squares, stripes and others [U EJ EJ [7]. Those works typically have 
been based on theoretical adsorption models calculated for convex particles like spheres, 
spheroids [8] and spherocylinders [9]. However, the majority of protein molecules are not 
convex and, therefore, there has been recent interest in the packing of concave objects 
[TU| inj [T2| [T3| [TS] . It appears that particle shape and symmetry can be crucial for 
fundamental properties of adsorbed monolayer [15], [161 El HE]- In this paper we study 
the irreversible deposition process of a dimer - the simplest, concave particle that could 
model plenty of bio-molecules. 

The following section outlines the model, algorithm and parameters used during 
the simulation of the dimers deposition process. The next part contains results 
and discussion. It is focused mainly on the random maximal coverage ratio. The 
ratios estimation accounts for several theoretical models including blocking function 
analysis. This section also describes autocorrelations and orientational ordering inside 
the covering layer. The paper ends with a short summary. 

2. Model and simulation procedure 

The adsorption process typically takes place when colloidal particles diffuse close to the 
surface. Due to adhesion this process can create a film consisting of randomly adsorbed 
molecules. Here we are focused on irreversible adsorption producing monolayers of 
adsorbate. The most straightforward approach to numerically simulate these processes 
is molecular dynamics (MD). The advantages of MD are prediction accuracy and control 
over most environmental parameters like temperature and the diffusion constant. The 
main drawback is performance. For this reason we decided to use another method 
known as continuum Random Sequential Adsorption (RSA), which had been successfully 
applied to study colloidal systems [19]. It is based on independent, repeated attempts 
to add a dimer to film. The features of a single step of numerical procedure used here 
were the following: 

- a virtual particle was created and its position and orientation on a collector was chosen 

randomly, accordingly to the uniform probability distribution. 

- an overlapping test with adsorbed earlier nearest neighbours of a virtual particle was 

performed. This test bases on checking if a surface-to-surface distance between 
particles is greater than zero. 

- if there was no overlap the virtual particle was adsorbed and added to an existing 

covering layer. Its position did not change during further calculations, which 
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reflected irreversibility of the process. 
- if there was overlap the virtual particle was removed and abandoned. 

Attempts are repeated iteratively. Their number is typically expressed using 
a dimensionless time unit: 

to = (1) 

where N is a number of attempts, = 2irr 2 stands for the coverage given by a single 
dimer (figure [T]) and Sc is a collector area. It is worth noticing that there is at least 




Figure 1. Model of a single dimer. In the following considerations we treat r as 
a length unit (r=l). 

one RSA algorithm allowing faster saturation of the underlying collector [20]. It bases 
on tracing uncovered area and uses this information to decrease overlapping probability 
for forthcoming RSA attempts. However we did not decide to use it here, because it 
complicates analysis of the standard RSA kinetics, which we want to compare with 
previous works. 

In the case of our simulations, the adsorption was stopped after T = 10 5 to- Collector 
sizes used varied from 20 to 200 r. Simulations were performed using fixed boundary 
collectors as well as collectors with periodic boundary conditions. In the first case centres 
of both circles forming a dimer had to be inside a collector area. For each collector we 
get at least 100 covering layers. The coverage ratio 9 is calculated as follows: 

= n d ^, (2) 

where is a number of adsorbed dimers. Typical coverages for three different values of 
coverage ratio 9: are presented in figure |2j Looking at the above drawing for jamming 
coverage (9 ~ 0.55), one can notice that the dimers density seems to be higher near 
the collector edges. This is an undesirable effect as we are mainly interested in the 
coverage ratio for an infinite collector. In order to control the systematic error caused 
by finite size and boundary conditions, the whole process was simulated over different 
sized collectors with fixed boundaries as well as with periodic boundary conditions. 
Results shown in figure [3] suggest that the bias of a measure (deviation from a pure 
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Figure 2. Typical monolayer samples for three different coverages: 9 — 0.1, 9 = 0.3 
and 9 » 0.55. The collector side length was equal to 40r. Fixed boundary conditions 
were used. 



quadratic fit) is less than 1% for the largest collector we used, independently on the 
specific boundary conditions used. On the other hand, for small collectors with fixed 
boundaries, higher densities near collector edges can be successfully used for producing 
systems with interesting optical properties such as micro-lenses. 

Most of results discussed later in the paper were obtained using largest collector 
[L = 200) with fixed boundaries. We checked that use of periodic boundary conditions 
does not have any significant influence on presented conclusions. 




Figure 3. Number of adsorbed dimers versus collector side length. Signs represent 
values taken from simulations for fixed boundaries (triangles up) as well as for periodic 
boundary conditions (triangles down). Lines are quadratic fits: n(L) = 0.086L 2 + 
0.115-L + 2.59 (solid line for fixed boundaries) and n(L) = 0.086L 2 - 0.066L + 1.654 
(dashed line for periodic boundary conditions). 
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3. Results and discussion 

3.1. The maximal random coverage ratio 

The main purpose of this work was to determine the maximal random adsorption ratio 
for dimers and compare it with results obtained for hard circles (spheres) p2J EH EH E2] ■ 
That ratio should be specified for an infinite collector area and infinite adsorption time. 
Despite controlling the error due to finite collector size one have to deal with finite 
simulation times. Particularly, in the case of large collectors, it is not certain if there is 
any possibility of adsorption after the simulation time and therefore approximation of 
maximal coverage depends on the RSA kinetics model (see figJl]). There were plenty of 




Figure 4. Coverage versus dimensionless time (??). Dots come from simulation. Lines 
represent example kinetisc fits, solid: max — Bit) ~ t~~ x / 2 and dashed: max — Bit) ~ 
exp(— ct). 

previous works in this area [19j (2TJ [22j [23j [HJ [25j 12S] and the general conclusion is that 
asymptotically: 

e max - e{t) ~ r l ' D . (3) 

for irreversible deposition of circles or unoriented squares. D here denotes collector 
dimension. The situation changes when adsorbed particles are ordered J2TJ [22]. For 
example deposition of oriented squares for long enough time obeys the following relation: 

ft--*)- 2 ^- (4) 

In case of present work the planar orientation of particles were chosen randomly with 
an uniform probability distribution. However this symmetry could be broken because in 
close proximity of previously adsorbed particle there are more space for parallely aligned 
particles than for perpendicular ones. Therefore asymptotic ([3]) and (HI) were compared 
in fig[5j Although both relations approximates the experimental data well, the fit ([3]) 
is slightly better in terms of linear correlation coefficient. The values of 9 max can be 
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Figure 5. Coverage versus t~ x l 2 and \nt/t. Bold points are taken from simulation. 
Dotted lines represent linear fits: 8 rnax — 9 = 0.54748 — 1.3787 • i -1 / 2 and 9 max — = 
0.54493 - 16.512 • Int/t respectively. 

obtained by interpolation of t~ l l 2 and Int/t to 0. Here they are 9 max = 0.5475 and 
9max — 0.5449 respectively. The difference is located within a 1% systematic margin of 
error coming from finite collector size. 

In order to compare our result with experiments, we analysed data obtained for 
adsorption of insulin. Assuming insulin dimer has a mass of 11616[-Da] and a size of 
875[A 2 ] one can find that 8 max = 0.55 corresponds to a monolayer surface density of 
l.2l[mg/m 2 ]. Typical values from experiments are in the range of 1.3 — l.Q[mg/m 2 ] 
depending on concentration and insulin type (human, Zn-Free) [271 EH] • Difference is 
noticeable but not significant. It can be explained by two causes: a more regular shape 
of the insulin dimer compared to our model (figured]) and the fact that, depending on 
concentration, the insulin particle can also appears as a monomer or hexamer. Both of 
them follow to higher coverages. 
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3.2. Blocking function and fluctuations in the number of adsorbed particles 

In the real experiment, adsorption kinetics depends typically on two factors: efficiency 
of the transport process (mainly diffusion or convection depending on experimental 
setup) that brings adsorbate from the bulk to the surface and the probability of catching 
particles, which are in a close proximity [El E31 [251 ESI ESI EQj EH [32], [33] . Here we would 
like to focus on the second factor, which is described by the blocking function, also 
known as the Available Surface Function (ASF). It can be easily obtained from the 
simulation as a ratio of successful attempts to all RSA attempts. Obtained ASF is 




0.1 0.2 0.3 0.4 0.5 

e 



Figure 6. The blocking function - successful attempts ratio versus coverage. Dots are 
simulation data, whilst solid and dashed lines are fits: ASF(6) = 1 — 4.770 + 5.586> 2 , 
and ASF{9) = (1 — 8/0. 543) 2 ' 97 correspondingly. The inset presents the same data in 
a logarithmic scale. 



presented in figure El The quadratic fit is commonly used to estimate adsorption rates 
at a limit of small coverages: 

ASF{6) = 1 - d9 + C 2 6 2 (5) 

In the case of dimers, simulations show that C\ = 4.77 and C-i = 5.58, whereas those 
parameters for hard circles adsorption are Cx — 4 and C 2 = 3.31 [25J. It shows that 
available surface shrinks faster for dimers. This follows intuition, because for successful 
adsorption a particle should have not only appropriate coordinates but also the right 
orientation. 

Moreover, ASF provides another way to estimate maximal coverage by analysing 
an adsorption probability. Specifically, the second fit in figure [6] suggests that maximal 
coverage is equal to 6 max = 0.543 and also provide additional support for model (j3J) 
because: 

<"«» - Asnm - d - I^V (6) 



dt x w/ V 9 



max , 
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yields 

e max - o{t) ~ r 1 ' 2 (7) 

Although estimating adsorption kinetics through ASF is straightforward and 
commonly used, there is a problem with direct measure of the blocking function during 
experimentation. Therefore, researchers are likely to count density fluctuations in the 
number of adsorbed particles for a given coverage. Typically, those fluctuations are 
expressed in term of reduced variance of particle number n inside a specified area: 
a 2 = o~ 2 (n)/ (n). It can be shown that, at least in the limit of small coverages, 
& 2 (9) = ASF {9) [31]. This comparison, in the case of our simulations, is presented 
in figure As expected, cr 2 (9) follows ASF [9) only for small coverages. For 9 > 0.2 

0.8 
0.6 

el 

ID 0.4 
0.2 

°0 0.1 0.2 0.3 0.4 0.5 

e 

Figure 7. Reduced variance of adsorbed particles as a function of a coverage. Dots 
correspond to the simulation data, whilst solid line is a quadratic fit J5]) to ASF. 

the difference between them becomes significant. 
3.3. Autocorrelations 

Autocorrelation of particles is another important characteristic of the monolayer. Here 
we are able to measure at least two different autocorrelation. The first is a standard 
distribution of distance between two molecules. Results are shown in figure El Lets 
first concentrate on reference frame given by autocorrelation inside monolayer build of 
spherical particles (pluses in figure [8]). The function has a maximum for r = 2.0 (nearest 
possible neighbour), then due to excluded volume it approaches minimum. Next, a much 
weaker maximum is around r = 4. Because of the random character of coverage, those 
oscillations vanish superexponentially [20J and functions stabilises at a value of 1, which 
is a result of normalisation. 

Autocorrelation for dimers (circles) is different. Function rises very slowly with 
distance to approach its first maximum around r = 3.8. It suggests for example that 
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Figure 8. Autocorrelation functions. The data were analysed in two different 
ways. Dots represents results for dimers as themselves whereas triangles shows 
autocorrelation for circles, regardless of that they are forming a dimer or not. Pluses 
show the reference results for circles. Inset shows asymptotic behaviour for small r. 
Fitted line is C(r) = -0.48151 • ln[(r - 2)/2] + 0.35683. 



there are not many dimers lying side-by- side, at a distance close to r w 2. Then, 
correlation approaches wide and shallow minimum around r = 5. This behaviour is 
mainly an effect of a dimers shape. We expect that for more anisotropic molecules, for 
example fibrinogens, there could be no minimum at all. 

Although at first it seems that autocorrelations for circles and dimers are totally 
different they can be easily compared when dimers coverage is treated as it was build 
up of independent circles. Here (triangles) autocorrelation looks almost the same as 
in the case for circles, only that maxima and minima are sharper. This similarity 
could explain why maximal random coverage for circles and dimers are so close. The 
asymptotic behaviour is also similar. At small distances it has the universal form derived 
in [2HE2!: 

C(r) ln[(r/2-l)] for r -> 2 + , (8) 

whereas for large r the decay seems to be even faster than for circles. 

3.4- Ordering 

The non-uniform shape of a dimer gives the possibility to check if any orientational 
order appears in a monolayer. Such ordering was widely investigated before, but mainly 
using lattice topology of collector surface eg. [35]. As mentioned in sec. 13.11 it could 
also influence on kinetics of RSA. 

To measure orientational order in our, continuous system we introduce the following 
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function determined by a dimers configuration: 

1 N 

S ( ( f ) ) = T^H( a; i cos + 2/i sin 0) 2 > ( 9 ) 

1=1 

where (xi,yi) are coordinates of a unit vector along the z-th molecule in a layer. It 
can be noticed that S(4>) is an average scalar product between molecules orientation 
and the direction given by an angle <fi. Thus, for an ideally aligned layer, S(<f>) will 
oscillate between and 1, where the maximum corresponds to an angle being parallel to 
molecules and minimums are reached for perpendicular to the direction of alignment. 
For pure random alignment S(4>) will be constant and equal to 0.5. In general, the mean 
orientation given by maximum S(<f)) can be estimated from: 

J2 N x ' 

taia(f) ex = — ^ 1=1 _ l ' . (10) 

The above equation is fulfilled by both maximum <fi max and minimum 4> m i n = <p max + 7r/2. 
Dependence between maximal value of S(<p) and collector size is shown in figure |9j 
Global order, although not very high, are strongest in small collector coverages. There 
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Figure 9. Order S(4> max ) dependence on collector side size. Trangles represent 
simulation data for fixed boundaries (triangles up) and periodic boundary conditions 
(triangles down). Dashed line is a simple analytical fit: S(<fi max ) = 1.58/L + 2.45/L 2 + 
0.5. 

could be at least two reasons for this. Firstly, if the allowed space for placing the 
following dimer is restricted, it is more probable to find enough room for parallel 
alignment than for a perpendicular one. Secondly, our adsorption conditions prefer 
parallel alignment at collector borders, because the dimer is placed down only when 
centres of two of their circles are touching the surface. For small collectors both, 
described above, effects are stronger however influence of fixed boundaries is irrelevant 
(figEJ). In order to determine which one of them is more important, we analysed local 
ordering in large collectors. It was done by calculating the mean value of scalar product 



RSA Study of Dimers 



11 



0.54 



A 




O 
V 




0.5- 







2 



4 

r 



6 



8 



Figure 



10. Correlation of ordering versus distance for collector size 200x200. 



between two dimers at a given distance. Results are drawn in figure [TDJ The local 
ordering and its range are quite small. It practically disappears when distance exceeds 
0.5. It suggests that ordering visible in tiny collectors is mainly due to their small size. 

3. 5. RSA of dimers in higher dimensions 

The RSA of (hiper) spheres and shows some universal behaviour regardless of collector 
dimension [T6l 1201 I2T1 122] . The best known of them are Feder's law fl3]) and asymptotic 
relation for autocorrelation function (jSJ). Results described in previous sections shows 
that those relations are also valid for dimers in 2D. Moreover the maximal random 
coverage ratio agrees with value obtained for circles. It is interesting if this is only 
accidental coincidence or more general feature. To address this question we looked at 
RSA of dimers in 3D. As a collector we used cube having side length 35 with a periodic 
boundary conditions and dimer was modelled by two touching spheres. Presented results 
are obtained from 50 independent simulation runs. 

Figure [TT1 presents coverage kinetics versus t -1 / 3 . The kinetics ([3]) fits well to the 
data. The random coverage for dimers in 3D is 0.37 and within margin of error agrees 
with value 0.0381 obtained earlier for spheres [20J. 

Autocorrelation function is presented in figJT2]Again we observe the same behaviour 
as in 2D and as for spheres. The decay for larger distances seems to be superexponential. 
On the other side, for small r the asymptotic agrees with flSJ). The linear coefficient are 
approximately two times smaller than values published for spheres random adsorption 
both in 2D and 3D [20]. 

4. Conclusions 



The maximal random coverage for a dimers monolayer is 0.547 and within the margin of 
error is not distinguishable from results obtained previously for adsorption of spherical 
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Figure 11. Coverage versus t 1 / 3 . Bold points are taken from simulation. Dotted 
lines represent linear fits: 6 max - 9 = 0.36986 - 0.74547 • i~ 1/3 . 
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Figure 12. Autocorrelation function for spheres, regardless of that they are forming 
a dimer or not. Inset shows asymptotic behaviour for small r. Fitted line is 
C(r) = -0.51176 ■ ln[(r - 2)/2] + 0.5384. 



particles. Also monolayer density autocorrelations are almost the same. On the 
other hand, calculated surface blocking function (ASF) is significantly different (C\ 
and Gi coefficient). It suggests other kinetics of adsorption, however the Feder's law 
03]) is maintained. Density fluctuations can successfully estimate ASF only for small 
coverages. Orientational ordering is imperceptible for macroscopic collectors but could 
play significant role in micro scale. 

Presented results for both maximal random coverages and autocorrelation function 
suggests that random packing problems for dimers and for spheres is governed by the 
same rules for D > 2. 

This work was supported by grant MNiSW/0013/H03/2010/70. 
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